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INTRODUCTION  AND  SUMMARY 

'n-iis  final  report  documents  the  results  of  Contract  No.  N00014-76-C-0699 . 

Ilie  objective  of  this  study  was  to  determine  the  feasibility  of  constructing 
a nonlinear  lens  capable  of  operating  in  the  8 to  14  micron  region  of  the 
infrared  spectrum. 

Under  a previous  ONR  contract^  MCAIR  proved  the  feasibility  of  a unique 
nonlinear  lens  which  had  substantial  bandwidth  reduction  capability  for  use  in 
remote  viewing  systems.  This  lens  takes  advantage  of  the  "variable  acuity" 
characteristic  of  human  vision  to  reduce  the  amount  of  information  that  must  be 
transmitted  in  a wide  f ield-of-view  high  resolution  imaging  system.  A brief 
description  of  this  concept  is  presented  in  Appendix  A.  A laboratory  brassboard 
of  this  system  was  constructed  under  contract  to  the  Office  of  Naval  Research  i 

(ONR),  (2)  and  is  presently  under  final  testing  at  MCAIR. ^ | 

The  system  developed  for  ONR  operates  in  the  visual  wavelength  spectrum  and  i 

consequently  is  useable  only  during  daylight  hours.  The  goal  of  this  effort  was  to 
determine  the  feasibility  of  achieving  the  same  wide  field-of-view,  high  resolution 
capability  during  night-time  hours  by  operating  the  sensing  lens  in  the  8 to  14 
micron  infrared  wavelength  spectrum.  If  successful,  this  sensor  system  could  fully  ! 

support  human  vision  under  clear  night  conditions;  a capability  impossible  to 
achieve  with  conventional  FLIR  imaging  systems. 

A flow  diagram  of  the  overall  effort  is  shown  in  Figure  1.  It  consists  of  two 
phases,  an  analysis  phase  and  a design  phase.  In  the  analysis  phase  the  optimum 
parameters  and  constraints  are  determined  that  establish  requirements  for  the 
actual  lens  design  effort  which  will  be  accomplished  in  the  design  phase. 


Ft^ura  1.  Work  Flow  Diagram 
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A computer  controlled  optimization  analysis  indicated  t'..at  the  lens  should  be 
designed  as  an  F/2.0  with  an  optical  quality  represented  by  a point  spread  function 
in  the  focal  plane  which  has  a one  sigma  width  (radius)  of  16  microns.  This  lens 
produced  the  best  performance  within  established  size  constraints.  In  addition,  best 
performance  was  achieved  when  the  lens  was  coupled  to  a serial  type  scanner. 

A narcissus  analysis  shown  in  Figure  1 was  also  performed  to  establish  any 
possible  constraints  on  lens  design  due  to  internal  optical  surface  reflections. 

The  results  of  this  effort  indicated  that  the  worst  narcissus  returns  were  eight 
times  below  the  detector  noise  level  (NEP).  It  is  therefore  concluded  that 
narcissus  is  not  a problem  with  the  infrared  nonlinear  lens  configuration  selected. 

The  lens  design  phase  indicated  that  realization  of  a practical  configuration 
for  the  infrared  nonlinear  lens  is  feasible.  The  preliminary  design  has  an  8 inch 
clear  aperture  and  consists  of  only  4 elements.  It  is  very  similar  in  size  to  the 
visual  lens  1.  The  optical  performance  of  this  lens  exceeds  the  F/2.0,  16  micron 
blur  requirements  because  an  F/2.0  lens  with  a maximum  blur  sigma  of  12  microns  was 
achieved.  System  performance  predictions  show  that  it  is  possible  to  achieve 
.Minimum  Resoluable  Temperature  (MRT)  performance  of  better  than  0.05°  and  spatial 
resolution  as  high  as  0.95  cycles/mr  (1.8  minutes  of  arc)  which  is  adequate  tu 
support  human  dynamic  resolution.  All  element  surfaces  have  reasonable  curvatures 
and  no  problems  are  anticipated  in  fabrication. 


*Number  in  parenthesis  refers  to  the  corresponding  item  in  the  List  of  References. 
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ANALYSIS 

Tile  analysis  is  divided  into  two  parts;  the  paranetric  analysis  which 
develops  the  Size/Performance  data  necessary  for  tradeoffs,  optimization,  and 
selection  of  design  parameters;  and  the  narcissus  analysis  which  develops 
constraints  for  the  optical  design.  These  efforts  are  discussed  separately 
in  the  following  sections. 

2.1  OPTiniZATION  ANALYSIS 

llie  overall  objectives  of  this  effort  are  to  establish  the  parameters  for 
an  infrared  foveal  lens  design  for  the  8-14  micron  spectral  region.  Since  specific 
mission  requirements  and/or  state-of-the-art  in  Infrared  variable  acuity  technology 
Is  not  defined  at  this  time,  the  approach  selected  was  a computerized  optimization. 
Tills  allowed  interactive  variation  of  key  parameters  while  size  and  performance 
were  assessed.  An  optimum  size/performance  balance  was  achieved  and  a preliminary 
lens  was  designed  for  these  optimum  parameters.  A flow  diagram  of  this  approach 
is  shown  in  Figure  2.  The  design  philosophy  is  as  follows;  Experience  from  our 
visual  lens  work  showed  that  the  lens  size  linear  dimensions  increase  inversely 
with  both  optical  quality  (blur)  and  F/number.  Conversely,  performance  as  measured 
by  Minimum  Resolvable  Temperature  (MRT)  improves  directly  with  these  two  parameters. 
Tlie  approach  shown  on  Figure  2 was  developed  to  determine  the  best  performance  but 
with  reasonable  size. 

Inputs  for  the  study  were  - 

(a)  Parameters  from  our  existing  visual  spectrum  lens 

(b)  Current  technology  scanner  and  detector  parameters 

(c)  System  inputs 

In  general,  these  Inputs  fall  within  two  categories,  those  that  remain  fixed 
throughout  the  study  and  those  tradeoff  parameters  that  will  be  tested  over  a 
range  of  values  to  determine  the  effect  of  each  value  on  performance  and  size. 

The  latter  are  F/number  and  blur.  A state-of-the-art  serial  scanner  and  a parallel 
scanner  were  also  evaluated  separately  in  the  analysis.  The  function  of  these 
inputs  are  described  below  along  with  rationale'  for  their  selection. 

LENS  OPTICAL  PARAMETERS 


Starting  at  the  upper  left  of  Figure  2,  the  optical  requirements  of  distortion, 
field-of-view  and  fabrication  constraints  are  used  to  define  focal  length  L(0)  and 
image  size  requirements  of  the  foveal  lens.  Basis  for  the  selection  of  the  value  of 
these  inputs  is  as  follows:  (1)  The  amount  of  distortion  used  was  made  equal  to 

the  existing  lens.  This  function  provides  a good  match  to  human  visual  acuity  and 
is  believed  to  be  near  the  maximum  distortion  technically  possible.  In  addition, 
any  sensing  system  constructed  with  this  function  will  be  compatible  with  the 
display  systems  constructed  on  related  contracts.  (2)  A field-of-view  goal  of  160° 
was  utilized  for  the  study.  This  value  was  selected  because  no  significant  reduction 
in  either  size  or  bandwidth  is  achieved  for  smaller  90°  - 100°  field-of-view.  In 


r ™ 


addition,  this  FOV  does  not  impose  any  difficulty  in  design.  Fabrication  of 
our  existing  visual  spectrum  lens  has  shown  that  the  rear  element  spherics  are 
extremely  difficult  to  produce  if  the  diameter  is  less  than  0.75  inch. 
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Figure  2.  Optimization  Analysis  Technical  Approach 
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The  relationships  established  by  these  parameters  are  shown  in  Figure  3. 
These  curves  were  used  to  select  an  on-axis  focal  length  of  2 inches.  This  then 
established  an  image  size  of  0.72  inches  for  160°  coverage.  A curve  for  a 100° 
FOV  is  also  shown  on  Figure  3 to  emphasize  the  point  made  above  that  little  is 
gained  through  FOV  reduction. 
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Image  Height  • In. 


Figures.  Focal  Length/Image  Size  Parametrics 


OPTICAL  RELAY  PARAMETERS 

I'/hen  the  image  size  is  known,  it  is  possible  to  define  the  optical  relay  that 
is  used  to  couple  the  image  to  the  detector  plane.  The  optical  relay  must  be 
defined  because  it  establishes  how  the  nonlinear  lens  F/number  appears  to  the  detector 
which  is  required  for  sensitivity  calculations,  and  to  establish  how  the  detector 
dimensions  translate  to  the  NLL  focal  plane  which  is  required  for  spatial  frequency 
calculations . 

Review  of  available  manufacturer's  literature  on  the  serial  and  parallel 
scanners,  makes  it  apparent  that  the  serial  type  placer,  the  most  severe  constraints 
on  NLL  design.  This  configuration  is  shown  in  Figure  4.  The  constraint  is  imposed 
by  the  high  speed  scan  drum.  This  necessitates  a stationary  exit  pupil  at  the  scanner 
facet  location.  The  optical  equivalent  of  this  type  scanner  is  also  shown  in 
Figure  4.  The  relationships  needed  for  analysis  are  computed  as  follows.  The 
collimator  focal  length  (L^)  must  be  sufficient  to  make  the  scan  cover  the  image 

which  has  a half-height  (H).  Then  the  angular  coverage  (6^)  is 


-£  = IL 

2 L, 


L = 

C tan  6 


1 


Tlie  detector  focal  length  (L^)  is  chosen  to  give  the  correct  resolution  element 
size  at  the  detector  (a  ) of  an  element  in  the  image  plane  of  size  (Ah), 


(3) 


Solving  for 


(4) 


Optical 

Equivalent 


OP77  032I-1 

Figure  4.  Optical  Relay  Parameters 
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The  detector  subtense  (Ah)  is  determined  by  the  number  of  scan  lines  (N)  required 

across  the  image  and  the  vertical  overscan  ratio  (q  ),  i.e. 

ovs  c 


Ah 


2H 

N ovsc 


Substituting  Equation  (5)  into  (4)  results  in 


a N L„ 

L = £ 

D 2 H n 


ovsc 


(5) 


(6) 


For  the  serial  and  parallel  scanners  under  consideration  the  parameters  are  as 
shown  in  Figure  5.  The  effective  F/number  (F#)  is  related  to  lens  F/number  by 


(7) 


This  is  the  effective  F/number  that  is  supplying  optical  energy  to  the  detector. 
It  will  be  used  for  performance  calculations. 


Type 

(in.) 

ay 

(in.) 

N 

'Iqvsc 

be 

(in.) 

•-D 

(in.) 

Serial 

0.0013 

0.0015 

488 

1.5 

1.29 

1.0 

Parallel 

0.002 

0.002 

360 

1.0 

2.46 

2.67 

QP77-Q32t-24 


Figure  5.  Scanner  Parameters 


Other  parameters  related  to  the  scanner  and  detectors  are  scan  efficiency, 
the  detector  geometry  and  detectivity  and  are  listed  in  Figure  6.  The  detectivity 
values  include  improvements  due  to  effective  cold  shielding.  Manufacturer's  actual 
data  was  used  to  estimate  the  Improvement  Instead  of  the  theoretical  values. 
Detector  equivalent  electronic  response  BU  was  determined  from  manufacturer's 
literature.  This  parameter  was  found  to  be  significant  only  in  the  case  of  the 
serial  scanner. 
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Serial 

Parallel 

Vertical  Scan  Efficiency 

93% 

N/A 

Horizontal  Scan  Efficiency 

82.7% 

75% 

Number  Parallel  Detectors 

] 

180 

Number  Serial  Detectors 

9 

1 

Detectivity  - D* 

5 X 10’° 

5 X 10’° 

Detector  Equivalent  Electronic  Response  BW 

13  (106) 

N/A 

OP77-032t-2S 


Figure  6.  Other  Scanner  Parameters 


LENS  SIZE  PARAMETRICS 


During  initial  phases  of  the  study,  the  lens  size  parametrics,  i.e.,  size 
as  a function  of  F/number  and  optical  blur  of  the  existing  visual  spectrun' 
lens  was  used.  These  parametrics  are  shown  in  Figure  7.  As  actual  lens  design 
progressed,  these  parametrics  were  updated.  However,  this  was  not  found  to  be 
necessary.  All  optimized  designs  fit  these  data  very  accurately. 


Figure  7.  Size  Parametrics  for  a 2-In.  Focal  Length  Lens 


PERFORM^VNCE  CALCULATIONS 


System  performance  In  terms  of  MRT  was  computed  using  the  Night  Vision 
l.aboratory  (NVL)  FLIR  model  which  was  suitably  modified  for  the  remote  viewing 
system.  lliis  model  and  its  modification  is  described  in  Appendix  B.  Other 
than  our  tradeoff  parameters  of  F/number  and  blur  and  the  detector  and  scanner 
parameters  discussed  above,  this  model  required  input  of  the  parameters  listed 
below : 


(a)  Wavelength  Range  - 8.5  to  13  microns.  This  range  was  selected  for  its 
good  atmospheric  transmission  over  reasonable  path  lengths. 

(b)  Atmospheric  transmission  - Original  plans  were  to  use  0.5  but  were 
changed  to  1.0  to  reflect  system  performance  rather  than  the  extremely 
variable  mission  performance. 

(c)  Optical  transmission  - Based  on  state-of-the-art  techniques  a 0.5% 
reflectance  for  each  of  10  elements  (two  surfaces  each). 


(d)  (e)  System  electronic  response  - Since  the  NVL  model  includes  no  internal 
noise  sources  other  than  the  detector,  electrical  response  was  made 
insignificant  for  all  systems  by  assuming  flat  response  out  to  the 
highest  spatial  frequency  (f^)  of  interest  which  was  determined  by 

f > — rhere  is  the  instantaneous  field  of  view.  This 

^ FOV 

assumption  is  entirely  adequate  for  the  NVL  model  because  (1)  response 
is  limited  by  detector  geometry  and  (2)  noise  bandwidth  is  limited 
by  human  observer  factors  such  as  eye  integration  time  and  the  matched 
filter  function  (see  Appendix  B) . 

(f)  Display  blur  - TTie  characteristics  of  the  direct  view  display  were  used 
for  this  parameter.  This  unit  has  a blur  sigma  in  angular  terms  of  0.26 
minutes  arc.  This  is  about  10  times  better  than  the  anticipated  system 
resolution.  When  utilizing  the  NVL  model  it  was  most  convenient  to 
convert  this  to  linear  spatial  blur  at  the  same  focal  length  as  the 
foveal  lens  which  is  3.9  microns. 


(g)  Display  magnification  - The  remote  viewing  concept  inherently  has 
unity  magnification. 

(h)  Frame  rate  - A TV  compatible  frame  rate  of  30  f rames /second  is  utilized 
for  all  analyses  with  two  fields  per  frame  (2/1  interlace). 

(1)  Eye  integration  time  - the  NVL  recommended  value  of  0.2  sec  was  used 
for  this  parameter. 

(j)  Ter.-.perature  - The  usual  300°K  value  was  used  for  nominal  scene 
temperature . 

(k)  Emissivity  - An  emissivity  of  0.8  is  utilized  for  both  target  and 
background . 

(l)  Threshold  S/N  - The  NTL  value  of  2.25  was  used  for  this  parameter. 
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(m)  Vibration  sigma  - The  vibration  sigma  selected  represents  the  residual 
sensor  motion  after  stabilization.  A value  of  ly  was  used  and  this  does 
not  cause  any  significant  performance  loss. 

PARAMETRIC  ANALYSIS  AND  SELECTION  RATIONALE* 

Using  the  parameters  defined  above,  the  analysis  and  selection  proceeds  as 
follows.  For  each  scanner  type  the  tradeoff  parameters  are  Iterated  throughout 
their  ranges.  For  each  set  of  tradeoff  parameters,  MRT  vs  f^  and  lens  size  are 

printed.  Selection  then  proceeds  as  follows.  First  a specific  value  of  f^ 

is  chosen  which  is  the  spatial  frequency  at  which  the  optimization  will  be  made. 
TTien  the  data  is  screened  to  determine  the  remaining  tradeoff  parameters  (F/number 
and  blur)  that  yield  the  smallest  system  for  a range  of  MRT  performance  levels. 

This  procedure  is  shown  in  Figure  8.  The  spatial  frequency  is  then  varied  and  the 
entire  process  is  repeated.  These  data  were  used  to  locate  the  point  where  a 
reduction  of  MRT  is  not  warranted  by  the  corresponding  increase  in  size.  The 
knee  points  of  the  curves  represent  the  best  selection.  Beyond  these  points 
excessive  size  results  with  little  improvement  in  MRT.  Finally  the  MRT  curves 
for  the  selcted  systems  (one  optimum  for  each  design  f^)  are  plotted  as  shown 

on  the  lower  right  of  Figure  8.  A subjective  choice  is  then  made  as  to  which 
is  most  acceptable  from  a mission  performance  and  vehicle  installation  standpoint. 
The  results  of  the  analysis  will  now  be  presented. 


• Compute  MRT  (fy.  Fno'^ 


/ f - "n  Si?e  I 
!■  n " I 

I fx  WIRT  I 


• Plot  Minimum  Size 

System  for  Each  fx/MRT 
Combination.  Select  Best 
Lens  Parameters  Based  on 
ASize/AMRT 


• Plot  MRT  Curves 
for  Selected 
Optimum  Systems 


O»77-03M-14 


*ThiS  step  not  necessary  for  current  scanners  vvhere  N is  specifier! 

Figure  8.  Parametric  Analysis  and  Selection 
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RESULTS 

Size/Per formance  parametrlcs  of  the  IK  nonlinear  lens  are  shown  on  Figure  9 
for  spatial  frequencies  of  0,3,  0.5,  and  0.7  cycles/milliradlan.  The  optimum  systems 
for  each  of  these  frequencies  are  designated  on  the  curves.  Parameters  of  these 
"optimum"  systems  are  listed  in  Figure  10,  The  first  observation  from  these  curves 
is  the  advantage  of  using  the  serial  scanner.  Performance  is  significantly  better 
and  the  size  is  smaller  regardless  of  spatial  frequency.  For  this  reason  the  serial 
scanner  was  selected  for  the  lens  design  phase.  The  performance  of  the  three  systems 
coupled  with  a serial  scanner  Is  shown  in  Figure  11. 

It  should  be  noted  that  the  next  increase  of  spatial  frequency  to  1.0  cycles/mr 
results  in  a very  large  size  lens,  60  inches  or  larger.  This  leads  to  selection  of 
the  0.7  cycles/mr  system  with  a serial  scanner  as  the  best  choice.  Its  size  is 
about  the  largest  that  can  be  accommodated  in  the  manned  vehicle  application  that 
is  under  consideration  for  the  system.  Also  the  spatial  frequency  (0.7  c/mr)  is  verv 
compatible  with  human  visual  performance.  To  illustrate  this  human  dynamic  visual 
resolution  is  also  shown  on  Figure  11.  Note  that  only  the  0.7  cycle/mr  system  has 
performance  sufficient  to  support  this  requirement.  Based  on  the  observations, 
lens  parameters  of  F/2.0  and  16  micron  blur  were  chosen  for  the  design  chase  of  the 


— — Serial  Scanner 
““  Parallel  Scanner 


F/2  1 6 Systems 


MRT  Required  at  Spatial  Frequency  Shown  C 


Figure  9.  Scanner  Size/Performance  Parametrics 
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Figure  10.  Optimum  Systems 
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Figure  11.  Serial  Detector  Performance  as  a Function  of  Resolution 
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NARCISSUS  ANALYSIS 


Before  proceeding  to  the  design  phase  a discussion  of  narcissus  is  in  order. 
Two  types  of  narcissus  analysis  were  conducted  during  this  effort.  The  first  was 
conducted  to  devise  the  most  efficient  method  of  relaying  the  lens  image  to  the 
detector  focal  plane.  In  the  second  effort,  actual  narcissus  images  were  located 
and  their  effect  on  system  operation  was  predicted.  Details  of  these  analyses  are 
presented  below,  Basic  narcissus  theory  is  presented  in  Appendix  C. 

RELAY  NARCISSUS 


Concern  with  the  possible  requirement  for  a field  lens  at  the  nonlinear  lens 
focal  plane  necessitated  an  early  parametric  analysis  of  the  optical  relay  with 
respect  to  narcissus  effects.  Narcissus  is  caused  by  the  cold  detector  "seeing" 
Itself  as  a result  of  reflections  from  relay  optical  surfaces.  This  relay  is  shown 
in  general  form  in  Figure  12. 
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Tlie  tiinction  oi  this  relay  is  to  relmage  the  lens  image  of  tlie  real  world 
onto  the  detector.  This  would  be  an  easy  optical  task  if  it  were  not  for  the 
scanning  mechanism  which  essentially  translates  the  real  world  image  across  the 
detector.  ilie  scanning  mechanism  requires  that  all  ray  bundles  cross  at  a common 
plane  shown  as  the  Effective  Scanner  Stop  of  diameter  D on  Figure  12.  A problem 
arises  because  the  optical  output  of  the  nonlinear  lens  is  not  compatible  with  this 
requirement.  The  chief  rays  emerging  from  the  image  plane  are  not  "well  behaved" 
as  they  are  in  a conventional  lens  and  therefore  will  not  arrive  at  a common  plane 
after  collimation.  Tliis  necessitates  use  of  an  aspheric  field  lens  at  the  image 
plane.  ^^Jlllle  field  lenses  create  no  special  problem  in  visual  systems,  they  can  cause 
very  serious  "narcissus"  problems  In  the  infrared  spectrum.  The  geometry  of  narcissus 
is  also  shown  in  Figure  12. 


It  is  easy  to  see  how  optical  surfaces  near  any  image  plane  generate  the 
worst  possible  narcissus  effects  because  they  reflect  and  return  the  image  back 
through  the  system  with  maximum  efficiency.  The  field  lens  is  in  an  ideal  location 
to  create  this  proble.  For  this  reason  the  Narcissus  theory  developed  in  Appendix  C 
was  applied  to  this  relay.  For  this  study,  a reflecting  surface  was  theorized  with 
curvature  (R)  and  separation  (S.,)  from  the  nominal  image  location. 


Narcissus  theory  requires  that  the  location,  size  and  strength  of  the 
narcissus  image  be  k.nown.  This  requires  definition  of  the  distance  (S)  from  the 
narcissus  image  to  the  image  plane,  magnification  of  the  image  (M) , and  reflectance  (t) 
of  the  reflecting  surface,  as  shown  in  Figure  12.  The  first  tV70  can  be  determined 
from  geometric  optical  theory.  Using  the  definition  of  distances  shown  on  Figure 
12,  the  distance  (S)  and  magnification  (M)  are: 


S = 


where 


s 

2 2S,  + R 


(8) 

(9) 


and 


■-C  *2 


[L‘  + (Sj-Spa-L^-l^)] 


(10) 


I 

I 


The  narcissus  Image  power  on  an  on-axis  detector  of  area  (A^)  is: 
4 X Aj  X y 

Pjj  = —2 in  watts 

where 


(11) 


T “ Surface  reflectance  determined  from  optical  material  and  coating 
characteristics 
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DS  cdx 

X = (-5^ — or  K — — whichever  is  smaller 
y = or  M • whichever  is  smaller 

2Lj^  2 


cax,  cay  " Effective  cold  assembly  width,  heif^ht 


This  is  the  worst  case  value  which  will  occur  at  the  center  display  scan  position 
as  shown  in  Figure  12.  For  other  scan  positions,  narcissus  returns  are  much  more 
difficult  to  compute.  However,  the  maximum  scan  angle  beyond  which  the  narcissus 
return  vanishes  is  relatively  easy  to  predict.  This  occurs  at  the  angle  (6^ ) shown 

in  Figure  12  where  the  marginal  ray  from  the  detector  strikes  the  curved  reflecting 
surface  normally.  At  this  scan  angle  only  this  ray  returns.  At  greater  scan  angles 
the  reflection  does  not  pass  through  the  stop  and  thus  no  narcussus  exists. 

This  maximum  scan  angle  is 

D 

6,  = arctan  — =1 


or  image  position 


^ + R ^c  “ 

If  the  image  position  (x)  is  divided  by  the  format  width,  the  portion  of  the  display 
occupied  by  the  narcissus  image  can  be  established. 

The  above  theory  was  applied  to  both  scanners  postulated  in  this  study. 
Characteristics  of  these  scanners  and  their  required  relay  components  were 
computed  to  match  the  detector  plane  to  the  lens  focal  plane.  The  additional 
parameters  that  are  required  for  narcissus  calculations  are  also  listed  in  Figure 
13.  See  Appendix  C for  definition  of  these  terms.  The  probable  serial  scanner 
includes  a reflective  detector  surround  to  reduce  narcissus.  This  narcissus 
reduction  technique  is  also  employed  in  the  parallel  scanner. 
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Scanner  Type 

Effective  Cold 
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Dimensions  cm 

9 

NEP  watts 
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Length  ■ In. 
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(in.) 

(L) 

1 

®cax 

1 

^cay 

Aq  - cm 

Detector 

(Ld) 

1 

Worst  Case  Serial 

12.5x10'  ® 

8.05x10“" 

1.291 

1 

4 

1 

Probable  Serial 

0.116 

0.0066 

12.5x10“® 

8.05x10  " 

1.29 

1 

4 

D 

Parallel 

0.0087 

1.882 

2.5x10"® 

2.57x10“" 

2.46 

2.67 

4 

B 
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Figure  13.  Scanner  Parameters 


The  equations  developed  in  previous  paragraphs  were  applied  to  the  relay 
designed  for  each  scanner  type  to  determine  how  the  strength  and  the  size  of  the 
narcissus  return  varies  with  respect  to  the  image  location  (S^)  as  location  of 

the  field  lens  surface  and  its  radii  of  curvature  (R^O  are  changed.  Separation 

(S^)  was  varied  ± 2 cm  and  4 radii  were  used,  (±1  cm,  ±“) . The  results  are  pre- 
sented in  Figures  14  to  19.  Figures  ]4  and  15  show  narcissus  for  the  worst  case 
serial  scanner  with  positive  and  negative  field  lens  radii.  Note  on  Figure  14  that 
with  a flat  field  lens  (infinite  curvature,  = “) , separation  must  be  maintained 
greater  than  approximately  ±0.5  cm  to  make  narcissus  acceptable  (below  the  NEP) . 

If  the  field  lens  has  positive  radii,  the  negative  separation  must  be  increased. 

As  shown  in  Figure  15  negative  radii  on  the  field  lens  reverse  this  situation. 

The  upper  curves  in  Figure  14  show  the  raster  ratio  which  were  derived  for  an 

infinitesimal  detector  dimension.  As  such,  it  gives  an  indication  of  the  spreading 
or  expansion  of  the  detector  narcissus  image  with  detector  field  lens  separation, 
not  actual  narcissus  image  size. 

Ttie  data  shown  was  expected  except  for  the  negative  curvature,  negative 
separation  case  shown  in  Figure  15.  Here,  an  Inverted  image  occurs  at  the  detector 
plane  at  about  -1.8  cm  separation. 

Figures  16  and  17  show  the  narcissus  reduction  obtained  by  adding  a reflective 
surround  to  the  detector  assembly.  This  reflective  material  reflects  300°K  photons 
from  the  detector  mechanical  assembly  into  the  relay  optics,  thereby  making  the 
detector  cold  assembly  appear  smaller.  Theoretically,  this  technique  can  reduce 
the  size  of  cold  area  to  that  of  the  detector  material.  As  a practical  value  in  this 
case,  however,  about  seven  times  the  total  detector  array  area  is  assumed.  Results 
using  this  scanner,  show  narcissus  is  insignificant  for  any  configuration. 

Figures  18  and  ]9  show  narcissus  for  the  parallel  scanner.  Here  narcissus  is 
acceptable  except  for  separations  of  zero  or  one  cm. 

Since  a serial  scanner  was  selected  and  the  modified  detector  surround  is 
easily  added  if  not  already  part  of  the  scanner  design,  we  conclude  that  relay 
narcissus  is  not  a problem. 
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Figure  T6.  Serial  Scan  Assembly  with  Reflective  Detector  Surround  Narcissus  as  a 
Function  of  Object  Position  (S-|)  and  Rear  Surface  Curvature  (R) 
Negative  Case 
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Raster  Ratio 


Separation,  S-|  ■ cm 


Figure  17.  Serial  Scan  Assembly  with  Reflective  Detector  Surround  Narcissus  as  a 
Function  of  Object  Position  (S-|)  and  Rear  Surface  Curvature  (R) 
Negative  Case 
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Separation,  S-|  • cm 
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Figure  19.  Parallel  Scan  Assembly  with  Reflective  Detector  Surround  Narcissus 

Positive  Case 


NONUNKAR  LENS  NARCISSUS 


Narcissus  effects  due  to  elements  of  the  nonlinear  lens  were  calculated  by  the 
procedures  outline  in  Appendix  C.  In  order  to  accomplish  this  task  within  a reason- 
able cost  and  time  schedule  an  early  lens  design  was  utilized  whicli  is  sligiitlv 
different  than  the  final  design.  In  addition,  all  elements  are  considered 
splierical  with  radii  equal  to  their  on-axis  curvatures.  Neither  of  the  above 
approximations  are  belived  to  introduce  significant  error.  This  will  be  verified 
during  early  analysis  during  tlie  fabrication.  Results  of  the  narcissus  analysis 
are  presented  in  Figure  20  fot  each  scanner  type.  Tliey  were  obtained  by  scaling  the 
detector  forward  to  the  nonlinear  lens  focal  plane  using  the  relay  parameters  developed 
earlier.  These  data  show  that  acceptable  narcissus  (significantly  below  KF.P)  can 
be  obtained  for  either  serial  or  parallel  scanners. 

CONCLUSIONS 

The  analysis  indiciites  that  the  F/1.0  lens  with  16  micron  optical  blur  for  the 
8-14  micron  spectral  region  is  optimum.  This  lens  can  be  used  most  effectively 
with  a serial  scanner.  Narcissus  does  not  appear  to  be  a problem  if  used  with  a 
properly  designed  serial  scanner. 
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Figure  20.  Narcissus  from  Nonlinear  Lens  Surfaces 
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Section  3 


LKNS  DESIGN  AND  OPTIMISATION 


Ihe  lens  design  phase  proceeded  along  lines  shown  in  Figure  21.  The  effort 
started  with  a materials  survey  to  find  tlie  best  material  for  construction  of  tlie 
infrared  nonlinear  lens.  tiermanium  was  selected  for  tlie  severe  aspherics.  The 
remaining  elements  could  be  made  of  otlier  materials  based  on  dispersion  require- 
ments. The  primary  reason  for  the  aspheric  element  selection  was  advice  from 
the  lens  fabricator,  Frank  Cooke.  This  material  can  be  worked  as  easily  as  glass 
and  presents  no  fabrication  problems.  In  addition,  this  material  can  be  obtained 
with  very  good  index  liomogeneity  where  all  other  materials  are  questionable 
in  this  area.  Considering  the  thickness  of  the  large  aspheric  elements,  homogeneity 
is  of  upmost  importance.  Another  reason  for  selecting  germanium,  was  its  low 
dispersion.  Ihe  dispersion  is  nearly  an  order  of  magnitude  smaller  than  that 
of  any  other  applicable  material.  Low  dispersion  greatly  simplifies  optical 
design,  especially  in  an  unconventional  design  such  as  this  one.  Finally,  ger- 
manium has  good  stable  mechanical  and  thermal  characteristics.  The  only  dis- 
advantage of  this  material  is  its  high  index  of  refraction  which  does  not  impose 
any  restriction  on  design  but  does  require  sophisticated  optical  coatings  to 
maintain  high  transmission.  The  required  coating  technology  exists  and  does  not 
outweigh  advantages  of  employing  che  high  index  materials.  Wliile  the  design  began 
with  only  the  aspherics  made  of  germanium,  it  was  later  shown  that  because  of  its 
low  dispersion,  the  entire  lens  could  be  made  of  germanium. 

The  next  step  in  lens  design  was  to  initialize  the  existing  lens  computer 
program  with  parameters  of  the  IR  lens  and  establish  a first  order  design  for  a 
single  wavelength,  /in  index  of  4.003  was  used  with  an  8 inch  clear  aperture 
diameter  and  the  original  nonlinear  lens  distortion  function.  The  initial  design 
was  established  for  an  11  micron  wavelength  and  F/2.0.  The  target  in  the  next 
design  iteration,  was  a lens  with  an  optical  blur  of  nr  greater  than  16  microns  when 
imaging  over  tlie  entire  8-15  micron  spectral  region,  /ifter  an  acceptable  design  was 
achieved,  modifications  were  attempted  to  produce  a well  behaved  optical  output 
from  the  image  plane,  to  facilitate  optical  relaying  to  the  scanner  and  detector. 

In  general,  this  consisted  of  establishing  a fixed  exit  pupil  for  the  F/2.0 
ray  cone  output  by  using  aspheric  field  lens  surfaces.  The  design  was  iterated 
until  a fixed  pupil  was  achieved. 


/\fter  a successful  design  was  achieved,  it  was  noted  that  there  was  a variation 
in  exit  pupil  illumination  with  a change  in  field  angle.  This  indicated  that  a 
non-uniform  F/number  existed  ’ the  image  plane.  liHien  this  problem  was  examined 
it  was  discovered  that  F/number  was  very  difficult  to  define  on  the  nonlinear  lens 
because  of  the  disparity  that  exists  between  radial  and  tangential  focal  length. 

/\  detailed  radiometric  analysis  was  required  to  determine  how  to  define  tne  lens  F/ 
number.  Tliis  analysis  is  presented  in  Appendix  D.  Hie  result  showed  that  the 
best  way  to  define  F/number  was  to  define  entrance  pupil  geometry  and  compute  an 
equivalent  F/number  by: 


F/No  = I ^ 


r(9)h(6) 
ab  sin  9 


where  ■;  = 


a,  b = 


Object  field 
Focal  length 
Radial  image 
major,  minor 


angle 

in  the  radial  direction 
he  i gilt 

radii  of  the  entrance  pupil  ellipse. 
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Using  this  equation,  the  magnitude  of  the  F/number  problem  was  established.  It 
was  found  to  be  excessive  and  led  to  reiteration  of  the  entire  design. 

The  best  lens  design  achieved  during  this  contracted  effort  is  shown  in 
Figure  22.  Its  blur  and  F/number  performance  are  shown  on  Figure  23.  The  lens 
design  computer  data  is  included  in  Appendix  E.  Reiteration  and  further  optimiza- 
tion of  this  design  are  continuing  in  an  attempt  to  further  improve  F/number 
uniformity.  We  believe  it  will  be  acceptable  for  operating  with  the  selected 
serial  scanner  if  not  by  design  then  by  designing  for  a smaller  F/number  and 
stopping  down  within  the  relay  to  achieve  uniformity.  As  a last  resort  a radial 
filter  or  variable  thickness  element  coating  may  be  used  to  achieve  final 
uniformity.  However,  we  do  not  feel  this  will  be  necessary. 
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Figure  23.  IR  Nonlinear  Lens  Performance 


I 

\ 

I 


27 


it 


J 


SecLUm  4 


CONCLUSIONS 

Tho  analyses  and  design  presented  herein  conclusive! y show  that  an 
infrared  nonlinear  lens  is  feasible  and  practical.  Tlie  optimum  lens  for  use  with 
the  selected  current  generation  scanner  is  an  F/2.0  lens  with  an  8 inch  clear  aperture. 
A remote  viewing  system  using  tliis  lens  would  theoretically  have  tlie  performance 
shown  in  Figure  24.  Tlie  lens  was  designed  for  an  optical  blur  of  16  microns.  This 
is  a worst  case  blur  that  the  design  program  does  not  allow  to  be  exceeded  any- 
where in  the  image  plane.  At  many  points  on  tlie  image  plane  performance  is 
considerably  better  (See  Figure  23).  On-axis,  where  detailed  viewing  is  accomplished 
the  blur  Is  only  6 microns,  also  shown  on  Figure  24.  The  improvement  is  subs- 
tantial but  not  as  much  as  might  be  expected  by  examining  only  the  ratio  of  the 
blur  reduction.  The  reason  for  this  is  seen  on  Figure  23.  Here  the  MTF's  of 
each  element  of  the  system  are  shown.  Blur  is  only  one  of  the  many  degrading 
mechanisms.  Wliat  has  been  achieved  tiirough  the  optimization  process  is  the  best 
balance  between  these  MTF's.  Reducing  blur  without  changing  other  parameters 
leads  to  limited  Improvement  in  overall  performance. 

Other  conclusions  from  this  study  concern  narcissus  and  ease  of  fabrication. 
Basically,  naricussus  does  not  appear  to  be  a problem  if  a properly  designed 
scanner  is  utilized.  Fabrication  is  expected  to  be  free  of  unusual  difficulty.  All 
aspherics  are  relatively  mild  and  the  use  of  germanium  for  all  elements  permits 
fabrication  by  techniques  similar  to  those  used  in  our  visual  lens  fabrications. 
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Figure  24.  Serial  Detector  Performance  as  a Function  of  Resolution 
8 14  p Spectral  Region 
F '2  488  Effective  Scan  Lines 
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Response 


Figure  25.  Serial  Detector/Nonlinear  System  MTF's 
F/2.0  Blur  Circle  16^1  488  Line  8 14/j  Spectral  Region 
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Appendix  A 


BRIKF  DFSCRiPTlON  OF  'IllF  RFMOTK  VIEWING  SYSTEM  (RVS; 

The  RVS  concept  is  based  on  the  fact  that  the  human  visual  capability  can  be 
represented  by  a resolution  capability  of  about  130,000  elements,  provided  that 
these  elements  are  sized  non-linearly  according  to  the  acuity  function  as  sliown 
in  Figure  A-1.  An  image  with  this  characteristic  requires  only  about  2 MHz  video 
bandwidth  at  30  Hz  frame  rates.  In  comparison,  standard  techniques  would  require 
over  1,000  MHz  bandwidth  for  this  f ield-of-view  (180°)  and  resolution.  Even 
at  smaller  f ields-of-view,  the  bandwidth  saving  is  significant.  A comparison  of 
bandv/idth  requirements  for  varying  f ields-of-view  for  the  conventional  linear 
acuity  function  and  for  the  RVS  foveal  concept  is  shown  in  Figure  A-2 . Appro- 
ximately two  orders  of  magnitude  decrease  in  BW  is  achieved  with  the  foveal  system 
at  FOV's  greater  than  20  degrees.  In  order  to  mechanize  the  concept  described 
above,  a method  must  be  devised  to  generate  an  image  which  satisfies  the  optical 
requirements  of  the  eye.  The  RVS  concept  contains  a lens  system  that  creates 
optical  "distortion"  by  varying  the  spacing  of  the  angular  resolution  elements 
to  duplicate  the  acuity  function  shown  in  Figure  A-1.  This  process  is  illustrated 
in  Figure  A-3.  The  lens  transfer  characteristic  required  and  the  technique 
for  reconstructing  the  image  at  a remote  location  is  also  shown  on  this  figure. 
System  operation  is  as  follows: 

The  image  transmission  system  scans  the  photocathode  of  the  vidicon  or  photo- 
detectors of  an  imaging  array,  transmits  this  signal  to  the  remote  location,  and 
recreates  the  image  on  a CRT  or  light  valve  tube.  In  the  original  RVS  concept, 
the  distorted  image  is  expanded  using  a lens  system  with  a transfer  characteri- 
stic identical  to  the  sensor  lens  and  imaged  on  a spherical  screen  concentric 
with  the  nodal  point  of  the  lens. 

Obviously,  for  the  above  image  transmission  system  to  perform  adequately,  the 
optical  axes  of  both  the  sensor  and  projector  must  have  the  same  alignment  as 
the  viewer's  eye.  The  initial  RVS  system  concept  used  the  approach  outlined  in 
Figure  A-4.  The  position  of  the  projector  is  slaved  to  the  camera  by  a high 
accuracy  position  servo,  with  the  camera's  angular  position  commanding  the  pro- 
jector's position  relative  to  fixed  ground  station  reference  coordinates.  The 
viewer  at  the  ground  station  thus  has  the  same  angular  perspective  as  he  would 
if  he  were  located  in  the  remote  vehicle.  The  sensor  and  projector  must  also  be 
aligned  with  the  viewer's  foveal  axis.  In  the  original  concept  a Honeywell 
oculometer  was  employed  for  this  function.  The  oculometer  measures  the  angle 
between  the  eye's  foveal  axis  and  the  projector's  optical  axis.  This  error  signal 
is  transmitted  to  the  remote  vehicle  and  commands  the  camera  to  move  until  the 
angular  error  is  reduced  to  zero.  As  the  camera  moves,  the  projector  follows 
through  the  slaving  loop.  The  control  mode,  presently  under  study,  is  somewhat 
different,  however.  The  observer's  head  position  instead  of  his  eye  position  is 
utilized  to  point  the  remote  camera.  The  operational  difference  resulting  from 
this  simplification  is  that  when  the  viewer  uses  his  peripheral  vision,  he  must 
learn  to  rotate  his  head  towards  the  area  of  interest  rather  than  his  eyes.  A 
reticle  may  be  required  to  show  the  observer  the  location  of  the  highest  acuity 
area  of  the  display. 
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Figure  A-1.  Human  Eye  Characteristics 
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Figure  A-2.  Bandwidth  Requirements 


rrdMsrniltef 

Piocessor 


Vidicon 

180‘*  FOV  Bandwidth  at  30  Frames/sec 
Conventional  1000  MH/ 

Variable  Acuitv  ‘ 3 MHz 


Image  Tiansfei  Charactei istic 


Projection 

Lens 


GP76-1037  no 


Figure  A-3.  Electro-Optical  Schematic 
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Appendix  B 


/Vi'PLlCATLON  OF  THE  NIGHT  VISION  LABOR^VTOKY  (NVIJ 
■raEKMAL  VIEWING  SYSTEM  STATIC  PERFORMANCE  MODEL  TO 

TTiE  RVS 

It  was  suggested  that  the  NVL  Thermal  Viewing  System  Static  Performance  Model 
be  used  to  evaluate  the  performance  of  the  Remote  Viewing  System  (RVS).  However, 
repeated  attempts  to  convert  the  RVS  parameters  directly  to  the  NVL  model  have  led 
to  tlie  following  problem.  The  radial  distortion  function  of  the  foveal  lens  does 
not  lend  itself  to  an  IfTF  analysis  as  a function  of  object  field  angular  spatial 
frequency  as  called  for  in  the  NVL  model.  All  parameters  can  be  converted 
successfully  except  for  tlie  scan  velocity  term  because  a linear  raster  scan  on 
the  lens  image  plane  will  create  a variable  angular  velocity  and  variable  direction 
scan  in  the  object  field.  This  is  depicted  in  Figure  B-1.  Extreme  complexity 
results  when  attempts  are  made  to  convert  spatial  into  temporal  frequency.  This 
is  illustrated  by  the  rotation  of  the  f bar  pattern  in  the  lens  image  plane 
shown  in  Figure  B-1.  Given  enough  time,  an  analysis  could  be  made  in  a manner 
compatible  with  the  NVL  model.  However,  the  analysis  is  much  simpler  if  performed, 
not  in  object  field  angular  frequency  (cycles/milliradian)  but  in  spatial  frequency 
terms  (lines/millimeter).  For  our  purpose  of  optimizing  the  RVS  lens,  it  is 
simpler  to  work  in  terms  of  spatial  frequency  on  the  foveal  lens  focal  plane. 

This  simplicity  arises  because  seven  of  the  nine  Ml'F's  are  independent  of 
object  field  angle  at  this  foveal  lens  focal  plane  location,  and  the  scan 
velocity  is  undirectional  and  uniform  at  this  location,  thereby  making  easy 
conversions  from  spatial  to  temporal  parameters.  The  only  non-linear  conversions 
necessary  are  simple  geometrical  ones  which  translate  from  focal  plane  to  object 
field  and  display  space.  The  advantages  of  working  in  the  spatial  frequency 
terms  will  become  clear  as  the  analysis  is  developed.  In  the  following  develop- 
ment, the  NVL  model  approach  will  be  used  precisely  but  will  be  applied  in  the 
foveal  lens  focal  plane  as  a function  of  linear  spatial  frequency  (l/mm).  Parameters 
will  be  covered  in  the  same  order  as  they  are  in  the  NVL  Report  , which 
describes  the  model  in  detail. 

B. 1 MTF's 

1.  Optical  MTF  The  optical  MTF's  consist  of  a diffraction  MTF  and  a Gaussian 
MTF. 

(a)  Diffraction  In  angular  terms,  the  diffraction  MTF  is  referenced  as  Equations 
(9)  and  (10)  of  the  NVL  report: 


H (f  ,6)  = - Tcos"^  A - A(1  - A^)^'^^l 

opt  x - i- 


where 


A = AF^^f^/L(6) 


where  L(6)  is  the  equivalent  focal  length  which  changes  over  a 50/1  range  as 
object  field  angle  9 changes.  The  angle  6 is  the  absolute  angle  between  the 
point  of  interest  and  the  lens  optical  axis.  At  the  foveal  lens  image  plane 

f 


(B-3) 


3A 


Figure  B 1.  Scan  Distortion  Introduced  by  FOVEAL  Lens 


where  S is  the  image  plane  spatial  frequency  and  f is  its  object  field 
angular^equivalent  measured  along  the  scan  line  projection  in  the  object 
field  (_  direction  on  Figure  B-1).  Solving  for  f in  Equation (B-3) and 
substituting  tiiis  for  f^  in  Equation  (fi-2). 

A = X F,,  S (B-A 

A X 

Since  the  F/number  of  our  lens  is  constant,  the  diffraction  MTF  is  no  longer  a 
function  of  object  field  angle.  Thus  we  may  write  H which  indicates 

that  the  MTF  is  a function  of  the  independent  variab?i  S only.  Note,  however, 
that  conversion  to  object  field  angular  spatial  frequency  is  very  simple 
because  focal  length  is  constant  over  small  angular  increments  and  may  be 
determined  from 


f = S L(6)  (B-5 

X 

where  ..  is  along  the  scan  line  projection  in  the  object  field 
likewise 

f = S I.(e) 
w y 

where  w is  normal  to  the  scan  direction  in  the  object  field 


(B- 


(b)  “ A similar  simplicity  exists  here.  The  MTF  equation  with  the  angular 

terra  b ot  Kquation  (11)  oi  Re ference (B-1) replaced  with  its  equivalent  is; 


11  (f  ,0)  = exp 
blur  X 


- L(6)^  - 


The  toveal  lens  inherently  has  a constiint  spatial  blur  over  its  entire  focal  plane, 
so  that  the  sigma  (u)  of  Equation(B-7) is  a constant.  Substituting  Equation (B-5) into 
(B-7)  we  see  the  blur  MTF  simplifies  to 


(S  ) = exp  [-  2ir^  o^S 


blur'  X 


(B-8) 


Thus  this  MTF  like  the  diffraction  MTF,  is  no  longer  a function  of  object  field 
angle  because  the  focal  length  variable  has  been  removed. 

2.  Detection  fflF  - The  spatial  filter  MTF  of  the  detector  is  defined  as: 


SinCirf  Ax) 

iL  (f  ,0)  = V--  ^ — = Sinc(f  Ax) 

I)et  x’  TTf  Ax  ' X (B- 

X 

It  is  also  complex  in  our  system  because  the  angular  projection  of  the  detector 
into  the  object  field  (Ax)  in  this  equation  varies  with  absolute  object  field 
angle  (6),  Since  the  detector  height  is  still  uniform  at  the  lens  focal  plane, 
shown  in  Figure (B-2) as  Ah,  Equation(B-9) can  be  restated  as: 


Sin(iTS^Ah^) 
ttS  Ah 

X X 


(B-10) 


Again  the  MTF  becomes  independent  of  object  field  angle.  Note  from  Figure(B-2) that 
the  detector  height  (2ih^)  is  a function  of  detector  slze(a),  detector  system  focal 

length  (Ljj) , and  relay  focal  length  (L^,),  viz: 

Ah  = a (B-11) 

X X — 

There  are  two  ways  to  arrive  at  the  required  and  based  upon  whether  the  scanner 

or  detector  characteristics  are  known.  If  the  scanner  is  available,  horizontal  FOV 

(0  ) and  instantaneous  FOV  (a)  are  known.  Then  the  detector  focal  length  is 
s 


a a 

L = ^ 

D a a 
X y 


(B-12) 
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Location! 

Figure  B-2.  Scanner  Resolution  Element  Geometry 

On  the  forward  end  of  the  scanner  we  must  cover  the  foveal  lens  image  jjlane 

height  of  2h(9  ) which  subtends  the  angle  6 . Thus, 

° max  = 


h(0  ) 

max 


0 

tan  s 


(B-13) 


Solving  for 


h^e  ) 

max 

tan  6 /2 
s 


(B-14) 


If  the  detector  characteristics  are  known,  the  focal  lengths  are  a function 
of  detector  size  (4h)  projected  unto  the  image  plane  as  shown  in  Figure  (B-2). 
Detector  size  2.h  can  be  computed  directly  from  either  the  on-axis  resolution 
required,  the  number  of  scan  lines  required  across  the  vertical  FOv,  or  bandwidtli, 
response  restrictions  and  frame  rate  requirements.  The  focal  lengths,  and 

are  then  selected  to  make  the  detector  dimension  appear  as  the  required  Th  at 

the  foveal  lens  focal  plane.  The  scanner  horizontal  FOV  given  by  scanner  mechanics 

must  cover  the  total  image  plane  height  (2h(6^^^)).  In  either  case  the  detector 


37 


1 


>rri'  becomes; 


'Wt‘=x>  • ^ 


(11-15) 


Again  this  MTF  is  independent  of  object  field  angle. 

3.  Detector  Electronics  MTF  - It  is  in  the  MTF,  the  detector  electrical 
response,  that  we  get  into  real  trouble  trying  to  work  in  object  field 
angular  space.  For  a conventional  linear  optical  system,  a linear  detector 
scan  velocity  converts  into  a scaled  but  linear  angular  scan  in  the  object 
field.  This  is  not  true  in  our  system  as  was  shown  in  Figure  B-1.  A linear 
scan  in  the  x direction  on  the  image  plane  results  in  angular  velocities  in 
botii  0 and  6 directions  in  the  angular  object  field.  Both  of  these  angular 

X y 

components  are  nonlinear  functions  of  both  x and  y position  on  the  image 
plane.  Thus,  converting  from  spatial  frequency  to  temporal  frequency  becomes 
very  complex.  All  of  this  can  be  avoided  by  working  in  linear  spatial  plane 
terms.  If  the  scanner  has  an  angular  scan  velocity  £,  then  the  linear  motion 
of  the  instantaneous  FOV  on  the  foveal  lens  image  is 


V = 6L„ 
x C 


(B-16) 


The  conversion  to  temporal  frequency  (f;  is  therefore 


f = V S 

X X 

This  is  a constant  conversion  and  not  a function  of  time, 
electronic  MTF's  of  the  NVL  model  are  valid.  These  are 


(B-17) 


Therefore,  all 


»'Det(f) 


Elect 


HbCI) 

4.  Display  - The  RVS  display  is  the  inverse  of  the  foveal  lens,  which  results 
in  a conventional  linear  raster  generated  on  the  CRT.  The  CRT  has  a constant 
spot  size  and  the  expansion  optics  has  a constant  blur  at  the  object  focal  plane. 
Again  this  MTF,  if  derived  in  the  linear  spatial  plane,  will  not  be  a function 
of  object  angle.  If  the  optical  blur  and  CRT  spot  size  are  combined  and 
assumed  to  have  a Gaussian  MTF,  a composite  sigma  (c^)  results  and  the  MTF  is: 

where  r is  the  physical  ratio  of  format  sizes ; viz 

^ ^ \ens  image  (B-19) 

I SPLAY  CRT 
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By  contrasc,  if  this  were  accomplished  in  the  object  angular  plane,  the  MTF 
would  be  much  more  complex,  viz 

o 2,  ,2^  2-, 

- (B-20) 

0 T 

L(G)‘'M 

where  M is  any  system  angular  magnification  from  object  field  to  the  viewer. 

Aga  in  the  simplicity  is  obvious. 

5 & 6.  Stabilization  and  Eyeball  - The  remaining  two  MTF's  are  the  only 
two  that  are  not  simplified  by  working  In  linear  spatial  rather  than  angular 
terms.  First,  stabilization  tends  to  be  angular  input  to  the  system.  Using 
the  .'ITF  from  the  NVL  report; 


■4,os«x)  ■ “P 

Converting  to  the  foveal  lens  image  plane  results  in 

IL  (S  G)  - exp  [-PS  \(0)^1  (B-22) 

Los  X , X 


Similarly,  the  eye  views  the  display  in  angular  terms.  The  N'VL  MTF  is 


£quation(B-23)must  be  converted  to  the  foveal  lens  image  plane 


(B-23) 


(B-2A) 


In  conclusion,  seven  MTF's  have  been  simplified  at  the  expense  of  two  that 
have  been  made  slightly  more  complex  by  the  conversion  to  linear  spatial 
frequency. 


B. 2 NEAT 

NEAT  can  be  used  as  defined  in  the  NVL  document  since  it  is  independent  of 

spatial  frequency,  therefore  no  changes  are  required.  However,  it  is  a function  of  F/ 

number  and  a discussion  is  in  order  concerning  which  is  the  correct  F /number 

to  use  in  the  calculation.  It  is  not  simply  the  foveal  lens  F /number  or 

objective  lens  F /number  as  stated  in  the  NVL  document.  For  NEAT  calculations, 

the  F/number  must  represent  the  actual  ray  cone  supplying  the  detector. 

In  general  our  lens  will  be  limited  to  a fixed  image  height 

that  is  based  on  fabrication  considerations.  The  scanner  must  cover  this 
total  image  in  the  most  efficient  manner.  The  optical  geometry  of  the  lens 
and  scanner  can  be  reduced  to  the  basic  arrangement  shown  in  Figure (B-3). 

On  this  figure  the  ray  cone  that  supplies  the  detector  is  4^.  Note  that 

a maximum  ray  cone  ((^^)  that  the  detector  can  utilize  exists.  Obviously, 
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FOVEAL  LENS 


SCANNER 


the  best  design  would  have  <t’3=4’4-  This  is  not  always  possible,  however, 

because  the  lens  may  be  required  to  operate  with  available  scanners.  Therefore, 
it  will  be  necessary  to  define  both  of  these  ray  cones  and  utilize  the  smallest 
for  N’EAT  and  diffraction  MTF  calculations.  Before  processing,  note  that  ray 
cones  may  also  be  defined  by  F 'numbers  as  showai  in  Figure  (B-A). 
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Figure  B-4.  Ray  Cone  Parameters 
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4>  = 2 tan  ^(D/2L) 
and  F/number  may  be  defined  by 


(B-25) 


F 


# 


1 

2 tan  (|)/2 


(B-26) 


Therefore,  we  can  use  F/number  instead  of  $ to  define  ray  cones.  By  the  nature 
of  collimation,  the  ray  bundles  between  the  collimating  lens  and  the  detector 
will  be  of  constant  cross  section,  i.e.. 


And 


Therefore 


And 


^#2  ^#4 


^//4  ■ L°  ^#2 

F -^F 
^#3  ^#1 


(B-27) 


(B-28) 


The  larger  of  these  two  F/numbers  must  be  used  for  both  diffraction  and  NEAT 
calculations.  Now  note  that  is  the  F/number  of  the  foveal  lens  while  Fy^^ 

is  that  generated  by  the  scanner  internal  aperture  stop  (D^)  and  we  may  write: 


F 


//I 


(B-29) 


For  any  particular  scanner,  will  be  known.  From  Equation(B-28) the  smallest 

F/number  which  translates  into  the  largest  rav  cone  we  can  supply  the  detector 
with  is: 


which  is  also  obvious  from  inspection. 


(B-30) 
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For  our  ONR  work,  a scanner  has  been  specified  while  for  the  USAF  effort  only 
a detector  geometry  is  available.  The  following  derivations  will  try  to  accommodate 
both  arrangements.  If  the  scanner  is  specified,  and  will  be  known.  If 

is  not  given  it  can  be  calculated  from  the  scanner  Instantaneous  FOV  (a)  and 
detector  size  (a)  in  the  same  plane.  Then 


y 


(B-31) 


Also,  if  the  scanner  is  specified,  its  maximum  scan  angle  (0^)  shown  in  Figure 

(B-2).is  known.  Since  the  scanner  must  scan  the  entire  image,  6^  must  cover 

2h(9  ).  Then 

max 


h(6  ) 

max 

" 6 
tan  s 


(B-32) 


It  is  now  possible  to  compute  F^^^  from  Equation  (B-27). 


a tan  6 / . 

_y s_ 

Cl  h(0  ) 

y max 


(B-33) 


where  F..  is  the  foveal  lens  F/number.  Remember,  for  the  diffraction  MTF  and 
NEAT  calculations,  we  use  the  larger  o f ^//3- 

If  a scanner  is  not  defined,  one.  must  be  theorized.  This  can  easily  be 
accomplished  by  procedures  outlined  in  Lloyd  _ Generally,  scanning  hardware 

(mirror,  drum,  etc.)  are  selected  to  accomplish  a given  scan  angle  6^  (number  of 

facets).  Collimation  lens  focal  length  (L  ) can  then  be  computed  by  Equation  (B-32). 

L> 

From  the  required  scan  lines  across  the  maximum  image  height  ^(0^^^^),  the  apparent 

detector  size  (Ah)  at  the  image  can  be  calculated.  The  angular  instantaneous  FOV 
of  the  scanner  can  be  calculated  as  follows,  using  the  on-axis  geometry  of  Figure 
(B-2): 


(B-34) 


Ljj  can  then  be  calculated  by 


a 

L.  = ^ 


(B-35) 


can  then  be  defined  by  matching  the  lens  F/number 


"s  = F^^ 


(B-36) 


P 


B.  3 ^tRT  Calculatl oriH 


Tlie  follt  ing  MRT  equation  modifications  are  required  so  that  the  computation 
may  be  performed  in  linear  spatial  frequency  terms.  First,  in  the  SVL  MRT  equa 
must  be  replaced  by  the  apparent  detector  size  at  the  foveal  lens  image  plane, 
i.e.,  it  must  be  the  Ah  defined  on  Figure (B-2).  As  previously  demonstrated  in 
Equation  (B-11), 


Ah 

y 


(B- 


Also,  in  the  MRT  equation,  it  is  best  to  compute  the  Q integral  in  terms  of 
temporal  frequency.  This  eliminates  the  velocity  term  in  the  MRT  equation 
and  makes  the  Q integral  easier  to  compute.  The  Q integral  is  therefore 


o 

Of  these  terms,  only  H^,  the  transfer  function  for  a rectangular  bar  of  width  w 

has  not  been  defined.  This  transfer  function  is  in  linear  rather  than  angular 
dimensions,  i.e.. 


(B- 


where 


X 


(B- 


The  MRT  equation  written  to  show  the  dependency  of  two  variables  is 


MRT(S^, 0) 


SR'K-  NEAT 


Ah  s q(f,e) 

y, . X 

*^e  'OVSC 


1/2 

(B- 


This  equation  results  in  an  MRT  very  weakly  dependent  on  6.  To  obtain  the  MRT 
for  any  field  angle  0 we  convert  the  spatial  frequency  term  into  an  angular 

frequency  term  by  using  Equation (B-5)containing  the  focal  length  function: 


f = S L(9) 
b X 


Note  this  will  be  the  angular  spatial  frequency  in  the  scan  direction  (target 

bars  normal  to  the  scan  direction).  It  could  be  related  to  f and  f but 

X y 

this  does  not  appear  to  be  required  at  this  point. 
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Appendix  C 


i'llKOKY  ()[•■  NARCISSUS  ANAI.YSKS 


The  problem  of  tlie  IR  detector  "seeing"  cold  focal  plane  surfaces  which  are 
reflected  by  tlie  system  optics  in  addition  to  tlie  warm  target  is  labeled  the 
"Narcissus  Problem".  The  source  of  the  narcissus  problem  is  Illustrated  in  Figure 
C-1.  In  Figure  C-1  the  detector  assembly  is  shown  at  the  center  scan  position. 

All  areas  of  the  rear  clear  aperture  of  the  optics  are  flooded  with  photons  emitted 
from  the  300°K  structure  surrounding  the  detector  assembly.  Most  of  these  photons 
are  absorbed  by  the  lens  structure  or  pass  out  into  the  object  field.  A small 
number  however,  are  reflected  back  to  the  image  plane  by  the  various  lens  elements 
and  therefore  some  of  these  returning  photons  fall  on  the  detector  to  alter  its 
electrical  output.  Alternatively,  tlie  emitted  photon  flux  field  from  the  warm 
structure  may  be  thought  of  as  being  disrupted  by  the  cold  detector  assembly  which 
for  all  practical  purposes  has  no  emission.  Therefore  a cold  "hole"  is  introduced 
into  the  flux  field.  The  narcissus  problem  is  generated  when  this  hole  is  made 
to  translate  across  the  optical  image  as  illustrated  by  the  dashed  lines  in  Figure 
C-1.  This  occurs  during  the  mechanical  scanning.  This  translating  region  of  non- 
emission will  cause  a fluctuation  in  the  reflected  photon  distribution  at  the 
detector  plane.  If  this  spatial  fluctuation  is  smaller  than  the  total  scan,  an 
unwanted  signal  will  be  generated  during  the  scan.  The  generation  and  prediction 
of  the  narcissus  image  are  described  below. 

The  best  way  to  compute  significance  of  the  "narcissus"  energy  is  to  assume 
that  the  detector  cold  assembly  is  emitting  at  300°K  and  trace  the  reflections  that 
return  to  the  detector  plane.  The  spatial  distribution  of  reflected  portion 
of  this  energy  in  the  detector  plane  is  the  photon  energy  that  will  be  absent 
in  the  actual  situation,  hence,  a "cold"  output  will  be  obtained  from  the  detector. 
This  process  is  shown  in  Figure  C-2. 

At  the  "start  scan"  position  shown  in  Figure  C-2 (a),  three  points  of  emission 
on  the  detector  cold  assembly  are  shown.  In  the  most  severe  situation  (negative 
magnification)  the  reflected  energy  will  appear  on  the  opposite  end  of  the  scan. 

For  example,  if  the  optics  are  theorized  to  form  an  image  somewhat  in  front  of 
the  focal  plane  by  reflection,  the  three  ray  bundles  will  converge  forward  of  the 
focal  plane  and  expand  towards  the  focal  plane.  This  energy  forms  a "blur"  area 
in  the.  focal  plane  as  shown  in  Figure  C-2 (a).  Remember  this  area  represents  a 
region  of  "less  photons"  than  the  surrounding  area.  Since  the  detector  does  not 
lie  within  this  area,  its  output  is  not  influenced  by  it. 

In  Figure  C-2(b)  the  detector  cold  assembly  is  shown  at  1/4  scan.  Note 
that  here  the  "narcissus"  image  has  moved  towards  the  detector  but  still  has  not 
influenced  its  output.  In  Figure  C-2(c)  the  detector  is  shown  at  center  scan. 

Here  it  lies  directly  within  the  "narcissus"  area  and  will,  therefore,  have  less 
output  or  will  have  a cold  output  created  by  the  lower  photon  input.  As  the 
detector  completes  the  scan,  it  moves  outside  of  the  "narcissus"  area  as  shown 
in  Figure  C-2(d)  and  is  no  longer  affected  by  narcissus.  Tlie  total  result  then 
Is  a "cold"  image  resembling  the  narcissus  blur  area  at  the  center  of  the  FLIR 
display.  Tliis  is  Illustrated  in  Figure  C-3. 


To  compute  the  true  effect  of  narcissus,  a ray  trace  of  all  points  within 
the  detector  cold  assembly  should  be  conducted.  Sufficient  ray  tracing  must  be 
conducted  to  map  the  power  density  in  the  detector  plane.  Much  of  this  laborious 
procedure  can  be  avoided  by  a series  of  slmpllcations . For  example,  if  the  portion 
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For  ar^alysis  purposes  only.  The  Cold  Detector  does  not  emit.  The  background  would  emit  this 
energy  if  detector  assembly  was  absent. 

OP77  0J«7  4 

Figure  C-2.  Narcissus  Image  Generation 


5 


Figure  C-3.  Narcissus  Display  of  Figure  C-2 


of  the  narcissus  image  that  contributes  to  illumination  of  the  on-axis  detector  can 
be  established,  simple  optical  and  radiometric  theory  can  be  utilized  to  define  the 
maximum  strength  of  the  narcissus  effect.  A similar  procedure  can  be  used  to 
estimate  the  size  of  the  narcissus  image.  This  is  developed  in  the  following 
paragraphs . 


On  Figure  C-4(a),  emission  of  the  on-axis  point  of  the  cold  assembly  is 
illustrated.  If  narcissus  analysis  is  being  performed  on  a surface  forward  of 
the  scan  mechanism,  the  only  energy  that  will  be  reflected  back  to  the  detector 
is  confined  within  the  F/number  cone  of  detector.  This  cone  is  defined  by  the 
basic  optical  design.  From  elementary  geometric  optics,  the  energy  returning 
from  the  forward  surface  will  be  confined  to  this  cone  but  altered  by  the  optical 
magnification.  With  reference  to  Figure  C-4  it  is  easily  seen  that  the  blur 
circle  radius  is: 


R = S 


S 

2 M FNo 


(C-1) 


where  S = Distance  narcissus  image  is  in  front  of  detector  plane 
M = Magnification 

The  above  are  shown  in  Figure  C-4 (a)  and  (b)  and  is  the  narcissus  power  contribution 
for  a single  point  on  the  cold  emitting  assembly.  Obviously,  the  detector  will  be 
affected  by  the  integrated  results  of  the  emission  of  all  points  on  the  cold  assembly 
that  can  pass  through  the  optical  system  and  be  reflected  back  to  the  detector. 

Three  of  these  are  shown  in  Figure  C-4(c).  It  is  easily  seen  that  all  of  these 
points  contribute  to  output  of  the  detector  located  near  the  optical  axis.  Tlie  difficult 
task  is  to  estimate  just  how  much  of  narcissus  image  is  capable  of  illuminating  the 
detectors.  This  is  a function  of  direction  of  the  chief  rays  forming  the  narcissus 
image . 


i 


In  Figure  C-5  the  image  cone  of  the  on-axls  detector  is  shown.  For  a correctly 
designed  optical  system,  it  will  be  impossible  for  any  returning  ray  to  reach  the 
on-axis  detector  unless  it  falls  within  this  cone.  Therefore,  rays  forming  the 
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Figure  4d. 
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Figure  C-4.  Narcissus  Blur  Circle  Radius 


N = 300°K  Radiance  of  Detector  Cold  Assembly 
T = Optical  Transmission  Reflectance  Product 
-N  = Radiance  of  Narcissus  Image 

IWq  = Detector  Area/S^  = Aq/S^ 

Pd=7NA,Wd 

a_  PD  = 4rNAoXY 
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Figure  C-5.  On  Axis  Narcissus  Power  Calculation 
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narcissus  image  cannot 
axis  than  the  distance 
is  therefore 

= 4 X y 


illuminate  this  detector  if  they  occur  farther  from  the  optical 
y on  Figure  C-5.  The  narcissus  area (A^)  on  Figure  C-5 


(C-2) 


Now  employing  another  simple  optical  theory,  that  the  image  radiance  (N^)  is  numer- 
icallv  equal  to  the  obiect  brightness  (N)  times  an  attenuation  factor  T,  which  is 
the  product  of  all  transmissions  and  reflectances  throughout  the  narcissus  path.  The 
derivation  of  the  on-axis  power  equaiton  is  shown  on  Figure  C-5.  For  a linear 
array,  the  last  equation  almost  always  applies,  i.e.: 


P 


D 


T N A^  M a 

T)  cax 


(C-3) 


Unfortunately  on-axis  power  is  not  sufficient  to  fully  assess  the  significance 
of  narcissus.  The  spatial  structure  on  the  display  must  be  known,  because  in 
many  cases  the  narcissus  power  may  be  many  times  the  detector  NEP  but  its  influence 
may  cover  all  detectors  equally  during  scan.  Therefore  it  will  not  be  seen.  This 
is  especially  possible  if  positive  magnification  exists. 


Starting  with  the  more  simple  case  of  negative  magnification,  the  case  depicted 
on  Figure  C-2,  estimation  of  narcissus  image  size  is  very  easy.  For  S > a , 
the  blur  image  width  (w)  in  the  detector  plane  is 


w 


2R 


S 

M FNo 


(C-4) 


Tlie  velocity  of  this  blur  during  scan  is 


V,_,  = MV. 

blur 


(C-5) 


Since  the  detector  moves  in  the  opposite  direction  with  velocity  V,  the  relative 
velocity  between  the  two  is 


V = V(1  - M) 
r 

The  time  required  to  scan  the  narcissus  image  is 

= 2R  ^ 2R 
'^blur  V(l-M) 


(C-6) 


(C-7) 


Hie  time  required  to  scan  the  total  format  is 


t 

s 


X 

V 


(C-8) 


where  X is  format  width 


The  ratio  of  the  display  width  covered  by  the  narcissus  image  is 

_ Slur  _ 2R 

t (l-M)x 

s 

For  positive  magnification,  tlie  theory  is  much  more  complex.  First  it  is 
necessary  to  know  the  nature  of  the  basic  image  to  proceed  with  this  calculation. 
For  a telecentric  image  (all  chief  rays  parallel  to  the  optical  axis).  The  tliaory 
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can  be  developed  as  follows.  Figure  C-6  must  be  used  for  this  derivation.  This 
figure  is  a view  in  the  x plane.  Therefore  the  detector  and  Narcissus  image  are 
vertical  lines  normal  to  the  paper  that  move  from  the  bottom  of  the  page  (start 
scan)  to  the  top  (end  scan).  In  this  case  (positive  magnification),  the  narcissus 
image  goes  in  the  same  direction  as  the  detector,  and  may  continue  to  illuminate 
all  detectors  throughout  the  scan,  unless  vignetting  occurs  within  the  optical 
system.  Again,  by  assuming  a telecentrlc  optical  image,  an  assessment  can  be 
made  of  this  effect  on  the  displayed  image.  In  this  case  the  optics  will  limit 
narcissus  rays  to  tlie  format  diagonal  dimension  on  the  detector  plane  because  the 
optics  were  designed  to  do  this.  This  means  that  no  ray  can  rise  higher  than  the 
height  (W)  in  the  narcissus  image  plane,  again  remembering  that  moving  vertical 
pencils  are  being  considered.  In  this  case,  the  narcissus  energy  will  be  completely 
vignetted  when  the  detector  array  is  at  a height  (— ) from  the  optical  axis.  It 
seems  logical  therefore,  to  assume  the  image  disappears  at  this  point.  The 
derivation  shows  that  for  this  assumption  the  ratio  of  narcissus  occupied  display 
format  height  is: 


r (S) 


M M a FNo 
cay 


(C-10) 


A similar  analysis  for  negative  values  of  S yields 


r (-S)  = 


M a FNo 
cay 
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Figure  C-6.  Narcissus  Size  for  Positive  Magnification 
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Appendix  D 

NON  LINEAR  LENS  F/NUMliER  DETERMINATION 


The  derivation  of  the  F/number  of  the  nonlinear  lens  is  defined  in  this 
Appendix.  Germaine  terms  are  defined  on  Figure  D-1.  Photons  are  radiated  from 
tlie  earth  (or  target)  at  a slant  distance  (3)  and  radial  distance  (P.)  from  the 
nonlinear  lens.  The  variables  S and  R are  much  greater  than  the  clear  aperture 
diameter  (DCA) . This  photon  flux  is  specified  by  a radiance  (N)  in  watts/CM^-Sterad ian . 
For  an  infinitesimal  area  in  the  object  field  (A  ),  the  power  (P^)  entering  the 
lens  pupil  of  area  (A)  is; 


where 


rhen 


P = NA  W 
a t 

uj  = Solid  angle  subtended  by  the  entrance  pupil 
A 


NA  A 
t 


(D-1) 

(D-2) 


(D-3) 


This  entire  flux  must  appear  on  the  lens  image  plane  within  the  image  (A^)  of 
object  area  (A^). 

A.  = dt  dr 
1 

The  power  density  of  the  focal  plane  is  therefore 


NA  A 
t 


A 9 

i S dt  dr 

Now  A^  can  be  related  to  the  image  plane  through  lens  focal  length.  By 
definition'" the  radial  focal  length  (f  ) is 

and  the  tangential  focal  length  (f^.)  is 


f J9) 


dt 

d9 


(D-4) 


(D-5) 


(D-6) 


where  6^  and  0^  are  the  radial  and  tangential  angular  subtense  of  A^.  Therefore, 

(D-7) 


dr  = f^(6)d(9^) 
dt  = f^(e)d(6|.) 

Now  A = Sd9  • Sd9  = S^d0  d0 
t r t r t 


(D-8) 

(D-9) 
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SubsLituLing  this  into  Equal  ion  (D-4)  we  have 


W 


N A d dh 
r t 


dt  dr 

Now  substituting  Equation  (h-8)  and  (D-9)  into  (D-10) 
N A 


W 


f^(e)f^(o) 


(IJ-IO) 


(D-ll) 


The  focal  length  fj.(9)  is  available  from  the  original  image  height  function  used 
in  lens  design  in  the  form  of  h(6)  vs  6,  thus 


f rni  dh(6) 


(D-12) 


The  tangential  focal  length  is  not  so  obvious.  It  must  be  calculated  by 
circular  symmetry  considerations.  The  same  number  of  area  elements  must, exist 
around  the  object  field  circumference  2itR  as  exist  around  the  image  plane  cir- 
cumference 2ir  h(9).  In  the  object  field  the  number  of  elements  is: 


^ 2ttR 

eo  d0  S 
t 


2tt  Sin 
d9 


(D-13) 


In  the  image  plane  the  number  of  elements  is 

2:rh(e) 


N . = 
el 


dt 


(D-14) 


Then  Equation  (D-13)  and  (D-14)  may  be  set  equal  to  each  other 


Sln9  _ h(9) 
d9^ 


dt 


(D-15) 


Or  by  definition  of  the  focal  length 
dt  _ h(9) 


= d9 


Sln9 


(D-16) 


Substituting  into  Equation  (D-ll) 
N A Sin9 


W = 


f/e)  h(9) 


(D-17) 


This  expression  can  be  evaluated  over  the  image  plane  as  a function  of  9 
to  determine  uniformity  of  focal  plane  Illumination  (W) . If  the  pupil  is  elliptical 
the  area  A is: 


A = •'Tab 


(D-18) 


Equation  (D-17)  then  becomes 
N~ab  Sin  9 


W = 


f^(0)h(0) 


(D-19) 
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To  put  this  in  relative  terms,  we  can  refer  to  the  equivalent  F/number  of  the 
nonlinear  lens,  i.e.,  the  circular  aperture  F/number  for  a conventional  lens  that 
yields  the  same  focal  plane  illumination  as  Equation  (D-19). 

If  A was  circular  with  diameter  (D)  and  focal  length  constant  (f)  the  flux 
through  A would  be 


where 


P = NA  — 


A = ttD 
4 


(D-20) 


^ NnP^  A^ 

4s2 


(D-21) 


If  the  lens  is  conventional,  all  this  power  falls  on  an  image  area 
f^ 

Therefore , 


i 4f  4(F  ) 

NOE-^ 

Equating  this  to  Equation  (D-19)  and  solving  for  F 


ff^(e)h(9) 

abSinS 


(D-22) 


This  equation  gives  the  nonlinear  lens  F/number  in  conventional  terms.  It 
is  used  to  assess  F/number  uniformity  of  the  nonlinear  lens  designs. 
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^ Appendix  E 

LENS  PREUMINARY  DESIGN  COMPUTER  DATA 

Titis  Appendix  contains  the  seventeenth  computer  run  out  of  50  runs  for  the  design  of 
tlie  nonlinear  lens  shown  scliemat  ical  ly  in  Figure  22. 
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